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Abstract 

Periodically structured materials can sustain both optical and mechanical excitations which are tailored 
by the geometry. Here we analyze the properties of dispersively coupled planar photonic and phononic crys- 
tals: optomechanical crystals. In particular, the properties of co-resonant optical and mechanical cavities 
in quasi- ID (patterned nanobeam) and quasi-2D (patterned membrane) geometries are studied. It is shown 
that the mechanical Q and optomechanical coupling in these structures can vary by many orders of magni- 
tude with modest changes in geometry. An intuitive picture is developed based upon a perturbation theory 
for shifting material boundaries that allows the optomechanical properties to be designed and optimized. 
Several designs are presented with mechanical frequency ~ 1-10 GHz, optical Q-factor Qo > 10^, mo- 
tional masses nietf 100 femtograms, optomechanical coupling length Lqm < 5 ^m, and a radiation-limited 
mechanical Q-f actor 2„, > 10^. 
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I. INTRODUCTION 



It has previously been shown that "defects" in a planar periodic dielectric structure can simul- 
taneously confine optical and mechanical resonances to sub-cubic-wavelength volumes^. As the 
co-localized resonances share the same lattice, and thus the same wavelength, the ratio of the op- 
tical to mechanical frequency of these modes is proportional to the ratio of their velocities. More 
recently, it was demonstrated that such co-localized resonances in a Silicon structure can strongly 
couple, via motion-induced phase modulation of the internal optical field, resulting in sensitive 
optical read-out and actuation of mechanical motion at GHz frequencies^. In this paper we aim to 
further develop the theory and design of these coupled photonic and phononic systems, laying the 
groundwork for what we term "optomechanical crystals". Here we choose a cavity-centric view- 
point of the interaction between photons and phonons, using the terminology and metrics from 
the field of cavity optomechanicJ^'^. An alternative viewpoint, more appropriate for guided- 
wave structures, may also be taken in which the interactions are described from a nonlinear optics 
(Raman-like scattering) perspective!^. 

We focus on two cavity devices in particular, a quasi-one-dimensional (quasi- ID) patterned 
nanobeam and a quasi-two-dimensional (quasi-2D) patterned nanomembrane, both of which have 
been studied extensively in the pasP^ for their photonic properties. The strength of the (lin- 
ear) optomechanical coupling in such structures is found to be extremely large, approaching a 
limit corresponding to the transfer of photon momentum to the mechanical system every optical 
cycleffSI. Simultaneously, the effective motional mas^ of the highly confined phonon modes is 
small, less than few hundred femtograms for a cavity system operating at a wavelength of 1.5 /jm 
and a mechanical frequency of 2 GHz. This combination of parameters makes possible the opti- 
cal transduction of high-frequency (multi-GHz) mechanical vibration j ^ l ^ l ^ •^ 1 ^ ^ I with near quantum- 
limited displacement sensitivitjff^^. Additionally, dynamical back-action^^ between the photon 
and phonon fields can be used to 

dampeiP™™ amphf^Pm 

mechanical motion, pro- 
viding an optical source of coherent phonon^25J26j y^,!^^^}^ ^^n then be used within other phononic 
circuit eiementP^imilll. Planar optomechanical crystals then, should enable a new generation 
of circuits where phonons and photons can be generated, routed, and made to interact, all on a 
common chip platform. 

Unlike the simple motion of a mirror on a spring in more conventional cavity optomechanical 
systemd^^I^, the complex mechanics of optomechanical crystal structures makes it difficult to in- 
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tuit the origin or strength of the optomechanical coupling. Nonetheless, understanding the nature 
of the coupling is crucial to the engineering of optomechanical crystal devices as the degree of 
coupling between different optical and mechanical mode pairs can vary by orders of magnitude 
within the same structure, with even subtle changes in the geometry inducing large changes in 
the optomechanical coupling. In the experimental demonstration of a nanobeam optomechanical 
crystal^, it was shown that the perturbation theory of Maxwell's equations with shifting material 
boundarie^ provides an accurate method of estimating the optomechanical coupling of these 
complex motions. Here we describe how this perturbation theory can be used to create an intu- 
itive, graphical picture of the optomechanical coupling of simultaneously localized optical and 
mechanical modes in periodic systems. 

The outline of the paper is as follows. We first analyze the quasi- ID nanobeam optomechanical 
crystal system. This nanobeam structure provides a simple example through which the salient fea- 
tures of optomechanical crystals can be understood. The mechanical Q of the structure is modeled 
using absorbing regions that provide a radiation condition for outgoing mechanical vibrations. The 
various types of mechanical losses are analyzed, and methodologies for minimizing or avoiding 
these losses are discussed. The dispersive coupling between the optical and mechanical modes is 
studied next. We use the aforementioned perturbation theory to analyze the optomechanical cou- 
pling strength, and which we display as an optomechanical coupling density on the surface of the 
structure. We use the density of optomechanical coupling picture to illustrate how the structure can 
be optimized to maximize the optomechanical coupling. Finally, we analyze the optomechanical 
coupling of a quasi-2D membrane structure, the well-known double-heterostructure photonic crys- 
tal cavit})^. We show how the optical and mechanical modes and their coupling can be understood 
in terms of the quasi-one-dimensional nanobeam example. 

II. ONE-DIMENSIONAL OPTOMECHANICAL CRYSTAL SYSTEMS: AN EXAMPLE 

To illustrate the nature of the optomechanical coupling and losses in OMCs, we will use a 
quasi- ID nanobeam structure which has been demonstrated experimentall)^^. Figure [T[ a) shows 
the general geometry of a periodic, quasi-one-dimensional OMC system made in a silicon beam 
of nanoscale cross-section. The system consists of an infinitely periodic array of by hy rectan- 
gular holes with center-to-center spacing. A, in a beam of width w and thickness, t (not shown). 
Although the actual structure will employ a defect to localize energy to a small portion of the 
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FIG. 1: (a) General geometry of the periodic nanobeam structure's projection (infinite structure, no defect), 
(b) Optical band diagram of the nanobeam's projection. The band from which all localized optical modes 
will be derived is shown in dark black, with Ey of the optical mode at the X point shown to the right of 
the diagram. The harmonic spatial potential created by the defect, along with the first three optical modes 
are shown as emanating from the X-point band-edge, (c) Mechanical band diagram of the nanobeam's 
projection. The three bands that form defect modes that will be discussed in this work are colored. The 
bottom-most mode is from the X point of the red band; the T points of the green and blue bands correspond 
to the middle and top mechanical modes, respectively. The frequencies of the defect modes that form from 
the band edges are shown as short, horizontal bars. 

beam, it is useful to consider the modes of this infinitely-periodic structure, since the structure has 
discrete translational invariance, allowing the optical and mechanical modes of the system to be 
classified according to their wavevector, kx, and a band index. We shall call the infinitely-periodic 
structure the projection of the system. The band picture provided by the projection allows a simple 
description of localized optical and mechanical modes as existing between two "mirrors" in which 
propagating modes at the frequency of the defect have a small or vanishing density of states; the 
mirrors surround a perturbation region where propagation at the modal frequency is allowed, local- 
izing the propagating mode between the mirrors. The optical and mechanical bands of the OMC's 
projection are shown in Fig. [TJb) and[ljc), respectively, for the structure A = 360 nm, w = 1400 
nm, hy = 990 nm, hx= 190 nm, and t = 220 nm. The material properties are parameterized by an 
isotropic Young's modulus, E = 169 GPa, and an index of refraction, n = 3.49. The optical bands 
are computed with the MIT Photonic Bands packageP^, while the mechanical bands are computed 
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with COMSOL Mutliphysicd^, a finite element method (FEM) solver. The structure does not 
possess a complete stop band for either the mechanics or the optics; nevertheless a defect in this 
structure can simultaneously produce highly-confined, low-loss optical and mechanical modes. 

The primary optical mode of interest will be the first TE-like (dominantly polarized in the y- 
direction) "valence" band mode at the edge of the first Brillioun zone (the edge of the first Brillioun 
zone is called X and the origin is called F). The electric field profile, Ey, is shown next to the band 
diagram. As described in detail in previous work on "zipper" optomechanical resonator^^^*^ using 
general momentum-space design rules of photonic crystal cavitied^, the localized modes that 
come from this band-edge mode are as far as possible from the light line while having a minimal 
amount momentum near = (and identically zero momentum at kx = 0) when used with a 
structure that is symmetric about a hole in the center. This reduces the radiation loss out of the 
structure. Because of the finite index contrast of the system, the optical Q is limited by radiation 
from optical momentum components that are close to k^ = 0, since the system can only guide 
momentum components that are above the critical angle for total internal reflection. This governs 
the design and choice of optical modes of the structure. Because the optical band has negative 
curvature at the X point, the frequency of the mode at the band edge must be increased to confine an 
optical mode coming from this band. This can be accomplished by decreasing A (making the holes 
closer together without changing the size of the hole). As has been showrP ^^ ^ '^'^l^^ l '^^ l '^^l^^ l ^^ l ^^ l 
both theoretically and experimentally, these nanobeam systems are capable of achieving very high 
radiation-limited Q-factors. 

Unlike light, mechanical energy cannot radiate into the vacuum. This makes the design rules 
for creating low-loss mechanical defect modes qualitatively different than those discussed above 
for optical defect modes, as will be discussed in the next section. Just as true photonic band- 
gaps are not necessary to achieve high confinement and low optical losses in nanowire structures, 
true phononic bandgaps are also unnecessary to achieve low mechanical losses. A quasi-stop- 
band, where a defect mode of a particular polarization, frequency, and fc-vector cannot couple 
a significant amount of energy to the waveguide modes of the mirror portion will be enough to 
acheive mechanical energy localization. Unlike in optics, all mechanical modes are "guided" by 
the structure, regardless of their ^-vector, which allows localized mechanical modes to be created 
from either of the high symmetry points, F or X. In fact, it will be shown that it is advantageous 
to draw mechanical modes from the F point, as this generally produces larger optomechanical 
coupling than drawing from X. Clearly the mechanical mode should be localized by the same 



5 



defect as the optical mode; if this is not the case, then the target mechanical band-edge should be 
essentially unaffected by the defect that creates the optical mode, and a separate defect must be 
found that can localize the mechanics without significantly affecting the optical mode. Finally, the 
localized defect mode that is formed from the band edge must be sufficiently optomechanically 
coupled to the localized optical mode(s) of interest. 
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FIG. 2: (a) Schematic illustration of actual nanobeam optomechanical crystal with defect and clamps at 
substrate, (b) Localized optical modes of the nanobeam OMC. The colors of the names correspond to the 
illustration of the inverted potential in Fig.[T|b). Localized, optomechanically-coupled mechanical modes of 
the nanobeam OMC. The colors of the names correspond to the colored bands and horizontal bars showing 
the modal frequencies in Fig.[TJc). 



The defect that will be employed to localize optical and mechanical energy to the center of 
the structure consists of a decrease in the lattice constant for the otherwise periodic array of Motai 
holes in the beam, as illustrated in Fig. ^c). For the modes that are localized by the defect, this 
effectively divides the structure into a "defect" portion where propagation is allowed, surrounded 
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by "mirrors", where the localized modes are evanescent, as discussed above. The particular defect 
used here consists of some odd number of holes, Nj^, with the spacing between the holes varying 
quadratically from the background lattice constant. A, to some value Ad, with the spacing varying 
symmetrically about the center hole (the hole dimensions are held fixed throughout the structure). 
The complete geometry, which we will refer to as "the nominal structure" is: Motai = 75, A = 360 
nm, w = 1400 nm, hy = 990 nm, = 190 nm, t = 220 nm, A^d = 15, and Ad = 0.85A. In 
the nanobeam structure described here, this defect simultaneously localizes many mechanical and 
optical modes. 

For localized modes, the quasi-harmonic spatial defect creates a quasi-harmonic potential for 
the optical and mechanical mode enveloped^. This creates a ladder of states for each band edge 
with approximately Hermite-Gauss spatial dependencies along the cavity axis (x), in direct analogy 
to the harmonic potential of ID quantum mechanics. As discussed above, the localized optical 
modes of interest come entirely from a single band-edge (the darkened band in Fig.[TJb)); the first 
three cavity modes of the defect from that band are shown in Fig.|2];b). Many localized mechanical 
modes with linear optomechanical coupling exist in this system. As examples, we will examine 
an exemplary optomechanically-coupled mechanical cavity mode from three different band-edges, 
even though each of these band edges produces a manifold of defect modes which may or may 
not have optomechanical coupling. The three modes are shown in Fig. |2jc), with the colors of the 
bands of Fig.[TJc) corresponding to the colors of the modes' label in the figure. 

III. MODAL CROSS-COUPLING AND MECHANICAL LOSSES 

Periodic structures can be fabricated to have phononic band gap^^SEHH^Iso]^ where mechan- 
ical energy loss by linear elastic coupling to the environment can be made arbitrarily small. 
Eventually, more fundamental losse^^^J^ such as thermoelastic los^SEMSIsl]^ non-equilibrium 
energy redistributiorP'^, phonon-phonon scattering^^, and the movement of dislocations and 
impuritied^^'^ should be accessible in these systems. First, however, the linear interaction of 
the optomechanical crystal and its surrounding substrate, which acts as a bath, must be understood 
and minimized. With this in mind, we use a finite element method model with weakly absorbing 
"bath" regions to model the losses in the system due to coupling of the mechanical energy into 
modes that are not confined. This method captures inter-modal coupling between the localized 
modes and all other mechanical modes of the system, some of which act as parasitic loss channels 
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FIG. 3: (a) In-phase and (b) in-quadrature mechanical displacement field (logjQ(|qp/max(|qp)) of the 
fundamental breathing mode of nanobeam OMC structure with weakly absorbing "pad", showing the prop- 
agating nature of the radiated mechanical waves in the pad region, (c) Dependence of Qm on the total length 
of the structure; the number of mirror holes on each side is {Nt — 15)/2. This shows the oscillatory of 
the pinch and breathing modes, which are coupled to waveguide modes, and the exponentially-increasing 
Qm of the accordion mode, (d) Mechanical band structure of the nanobeam OMC, with arrow tails indicat- 
ing the frequency and high-symmetry point of the breathing (blue) and pinch (red) modes, and arrow heads 
indicating the equi-frequency waveguide mode that acts as the dominant source of parasitic coupling. The 
effective bandgap of the accordion mode is shown in transparent green, with its frequency indicated as a 
horizontal green bar at the F point. 

into the surrounding "bath". 

The lack of a mechanical bandgap means that the superposition of fc- vectors necessary to create 
a localized mode in the defect coincide with ^-vectors of equi-frequency propagating modes of 
the phononic crystal mirrors. The localized modes and propagating modes of equal frequency will 



8 



hybrizide and couple whenever the symmetries of the modes do not forbid it. In addition to the 
propagating modes, there are "body modes" that exist purely because of the boundary conditions 
(and thus not represented in the band structure), such as vibrations, density waves, and torsions 
of the finite, clamped structure. If the simulated exterior boundary conditions allow energy in 
the propagating and body modes to be lost, the propagating and body modes that couple to the 
localized mode will act as parasitic loss channels for the localized mode. 

In a fabricated structure, the cantilever is attached to a substrate at both ends, rather than a 
hard boundary at the end of the cantilever. These more realistic boundaries must be included to 
model propagating and body mode losses. The propagating modes travel down the nanobeam and 
partially reflect at the contacts due to an effective impedance mismatch caused by the geometric 
change between the nanobeam and the bulk. The rest of the power radiates into the bulk, causing 
a loss of mechanical energy. Thus the locaUzed mode is coupled to a propagating mode with 
identical frequency that can radiate part or all of its energy into the surrounding "bath". This 
propagating mode also forms a coupled cavity resonance with the localized mode because of the 
reflections at the clamp points. The body modes have a softer boundary condition than q = at 
the boundaries, extending the body mode amplitude into the substrate. The part of the body mode 
that extends into the substrate can excite radiative modes of the substrate. The body mode then 
acts as a loss channel for any localized mode to which it is coupled. The localized, propagating, 
and body modes form a set of coupled resonators. Since the body and propagating modes are 
very sensitive to the total length of the structure, the self-consistent solution, which determines the 
loss of the localized mode, is very sensitive to the exact boundary conditions. Thus, to accurately 
simulate the true spatial profiles and losses of localized mechanical modes, one needs a simulation 
that reflects the true boundary conditions. 

To model the loss due to coupling to radiative modes of the substrate, we include a large, 
semi-circular "pad" on each side of the nanobeam, with the same material constants as the 
nanobeam. To make the pad act like a "bath", we introduce a phenomenological imaginary part 
of the speed of sound in the pad region; i.e., Vpad — vsiiicon(l + i^), where v = y^E/p. This 
creates an imaginary part of the frequency, and the mechanical Q can be found by the relation, 
Qm = Re{Vm} / (2Im{Vm}). By adding loss to the pad material, propagating modes will reflect 
part of their power at the contacts (the interface between the cantilever and the substrate) because 
of the change in the impedance from the absorption, not just the geometric change in impedance. 
From this point of view, T| should be made as small as possible, since this contribution to the re- 
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flection coefficient is an artifact of the simulation and is not present in the real system. However, T| 
must also be large enough that the self-consistent solution includes a propagating, radiated wave, 
which only happens if the wave is appreciably attenuated by the time it reflects from the edge of 
the simulation (where q = 0) and returns to the contact. Thus, the pad is made as large as possible, 
given computational constraints, and the absorption is increased until changes appreciably, 
which gives the threshold value for r\ at which the reflectivity of the contacts has an appreciable 
contribution from the absorption. The simulation is then performed with a value of r| that produces 
a propagating wave in the pad without causing an artificial reflectivity at the contact. Propagation 
in the pad is easily verified if the position of the nodes/antinodes swap between the in-phase and in- 
quadrature parts of the mechanical cycle (the nodes/antinodes of a standing wave are stationary). 
Figure ^a) and ^h) show the in-phase and in-quadrature (respectively) parts of the mechanical 
cycle of the breathing mode, clearly showing a propagating radiative mode in the weakly absorb- 
ing pad, with logio(|qP/niax(|qp) plotted to elucidate the attenuation of mechanical radiation in 
the pad. 

Limiting the artificial reflection at the interface of the non-absorbing and absorbing portions 
sets the maximum absorptivity, which in turn sets the minimum size of the pad (to guarantee 
that the radiation is completely attenuated before returning to the source). This size/absorptivity 
trade-off can be improved by making r\ vary as a function of position in the pad, starting out 
at zero and increasing radially outward (quadratically, say). This is analogous to a mechanical 
perfectly matched layer (PMLj^, which has the benefit of increasing the round-trip absorption 
while maintaining a minimum reflectance at the clamp due to absorption. This provides the same 
reflection-free absorption of mechanical radiation in a more compact simulation space, making 
better use of computational resources. 

Changing the length of the structure changes the resonance condition for both the propagating 
and body modes. This changes the amount of coupling to the localized mode in the self-consistent 
solution of the system. Thus the parasite losses into a waveguide mode should be periodic. Figure 
[Sj^c) shows Qm for the pinch, accordion, and breathing mode as the total number of holes (i.e., the 
length of the optomechanical crystal) is varied. For the breathing and pinch modes, the mechanical 
losses oscillate as a function of the total length of the nanobeam. For this particular geometry, the 
losses are dominated by propagating modes, and the oscillation period of the Q can be matched 
to a ^-vector of a waveguide mode in the band structure in Fig. [Sj^d) (shown as a dotted fine 
extending from the defect frequency). Thus the length of the structure can be tuned to minimize 
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mechanical losses in cases where a complete mechanical bandgap is not present. Interestingly, 
the Q of the accordion mode increases exponentially with the number of holes, indicating that 
the mode is evanescent in the mirror portions. Examining the four bands that cross through at 
the frequency of the accordion mode, we find that all four bands have a mirror symmetry about 
either the x — zoxx — y planes that forbids any hybridization or coupling between to the accordion 
mode. This creates an effective bandgap (shown in translucent green). Practically speaking, this is 
a much weaker stop-band than a true bandgap, because any defect in the structure that breaks the 
symmetry of the accordion mode about the mirror planes will cause a coupling to the waveguide 
modes in the gap. However, it is exactly this kind of symmetry-dependent effective bandgap that 
is responsible for the high optical ^2 | 38|45 | 46 J experimentally-fabricated structures. This 

gives some confidence that it is possible to fabricate structures that are defect-free to the degree 
necessary to achieve high Q^- 



IV. OPTOMECHANICAL COUPLING: DEFINITION AND INTEGRAL REPRESENTATION 

Cavity optomechanics involves the mutual coupling of two modes of a deformable structure: 
one optical and one mechanical. The optical mode is characterized by a resonant frequency (Oq = 
27rVo and electric field E(r). The mechanical mode is characterized by a resonant frequency = 
27rVni and displacement field Q(r) , where Q(r) is the vector displacement describing perpendicular 
displacements of the boundaries of volume elements. The cavity optomechanical interactions of 
the distributed structure and its spatially-dependent vector fields, E(r) and Q(r), can be reduced to 
a description of two scalar mode amplitudes and their associated mode volumes, with the coupling 
of the amplitudes parameterized by a single coupling coefficient, ^OM- 

The mode amplitude, c, and complex vector field profile, e(r), are defined such that the com- 
plex electric field is E(r) = ce(r) (the physical field is given by the real part of E(r)e"^0- For 
pedagogical reasons, the amplitude c is normalized such that the time averaged electromagnetic 
energy is equal to |cp; i.e. U = |cp = \ J dy£|Ep. This forces e to be normalized such that 

1 2 

1 = ^/ dVE|e| . In cavity quantum electrodynamics, one typically defines an effective optical 
mode volume, Vo = T dV ( — Y?^^ | , in order to gauge the strength of light-matter interac- 

\ max( V £E J / a a a a 

tions. 

The mechanical vibration's amplitude, a, and mode profile (displacement), q(r), are defined 
such that Q(r) = aq(r) . Here, a is defined as the largest displacement that occurs anywhere for the 
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mechanical field, Q(r), so that max(|Q(r)|) = 1. It is important to note that this particular choice 
of a determines the mechanical mode's effective volume and effective mass, Vm and meff = pVm> 
respectively. In order to represent an energy and be consistent with the equipartition theorem, this 
choice of a requires the complimentary definition meff = p/dv(jJ^) .To see this, note that 
the free evolution of the mechanical oscillator has, by definition, a time-independent total energy 
^^mechanical = ^(^^^cc^ + 06^). On the Other hand, integrating the total energy of each volume 
element must also give this same total energy. If we pick the point in phase space at which all the 
mechanical energy is potential energy (i.e. the classical "turn-around point"), we must have that 
^'mechanical = 5^^^/ PlQ(r)PdV = imeffH^a^, or, in other words, mefra^ = / p|Q(r)|2dV. One 
can arbitrarily choose the definition of the amplitude or the mass, but choosing one determines the 
other. Note that a is also the amplitude of zero-point motion of the canonical position operator 
in a quantized treatment. For a system like a localized mode of a phononic crystal defect cavity, 
where only a very small, localized portion of the total mass undergoes appreciable motion, the 
most sensible choice of the mass is the amplitude-squared weighted density integral, which, as 
stated above, is the choice of mass associated with a = max(|Q(r)|). 

The optomechanical coupling affects the optical mode by tuning the optical resonant frequency 
as a function of displacement, C0o((x); whereas the coupling affects the mechanical mode by apply- 
ing a force, which is expressed as a gradient of the cavity energy, d |c| /da. The optical resonant 
frequency is usually expanded in orders of the (small) displacement, a around some equilibrium 
displacement, ao. 



a)o(a) = a)o -F(a-ao)^ -^(a-ao) 

a=ao dec 



2 d^g)o 

a=ao ' ^"^ """^ da^ 



+ ... (1) 

a=ao 



In the case that the terms higher than first order can be neglected, this equation simplifies to 



tOo(a) = tOo 



+ (a-ao)— 
a=ao aCC 



where tOo = tOo 



= cOo + (a-ao)gOM = cOo+(a-ao)7^ , 

a=cxo Lqm 



is the equiUbrium resonance frequency of the optical mode, goM = ^ 



a=ao 



(2) 



a=ao 



is the derivative of the resonance frequency of the optical mode evaluated at equilibrium, and Lqm 
is the effective optomechanical length of the system. The effective length, Lqm, is a universal 
parameter that relates displacement to a change in optical frequency (i.e. a/ Lqm = ScOo/Wo)- 



From the definition, Loj;, = ^ 5^ 



a=ao 



• gou/ tt>o> one can see that reducing Lqm maximizes the 
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optomechanical coupling. It is simple to show thatLoM is equal to the spacing between the mirrors 
of a Fabry-Perot cavity when one mirror is allowed to move along the cavity axis or the radius of 
a microtoroid/microdisk for a radial breathing motion. For a "Zipper" cavity or double-microdisk, 
LoM is an exponentially decreasing function of the spacing between the coupled elements, with 
LoM approaching half a wavelength of light as the spacing approaches zero. 

The perturbation theory of Maxwell's equations with shifting material boundarie^ allows one 
to calculate the derivative of the resonant frequency of a structure's optical modes, with respect to 
some parameterization of a surface deformation perpendicular to the surface of the structure. If 
the result of a mechanical simulation is the displacement field, Q(r) = aq(r) = aQ(r)/max(|Q|), 
then 



where d = £e, n is the unit normal vector on the surface of the unperturbed cavity, Ae = £i — £2, 
A(e^^) = £]"^ — £2 ^ El is the dielectric constant of the structure, and £2 is the dielectric constant 
of the surrounding medium. 

To calculate Lqm by deforming the structure, one must simulate the fields with a deformation 
amplitude, a, that is large enough to be detectable numerically but small enough that higher or- 
der dispersion does not affect the frequency shift. To verify that higher order dispersion is not 
included, one must simulate the optical fields for a range of displacement amplitudes and extract 
the linear dispersion. Because perturbation theory can calculate the linear term exactly from a sin- 
gle calculation using the undeformed structure, this method has clear advantages over numerical 
methods using finite deformations. 

V. OPTOMECHANICAL COUPLING: VISUAL REPRESENTATION AND OPTIMIZATION 

In addition to being computationally simpler than deformation methods, the perturbative 
method of calculating the optomechanical coupling allows one to represent the optomechanical 
coupling as a density on the surface, with different parts of the structure contributing different 
amounts of optomechanical coupling. This yields much more information than just the value of 
Lqmj itself. The optomechanical coupling density is given by 




(3) 



ComW = -(q-n) [A£|ei||^-A(£ ^)\&^ 



2 



(4) 
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FIG. 4: For the fundamental breathing mode and the fundamental optical mode in the nominal structure, (a) 
FEM simulation of individual unit cell contributions to the total optomechanical coupling (each point com- 
puted by integrating ^om (Equation [3]| over the respective unit cell), (b) surface plot of the optomechanical 
coupling density, ^om- (c) surface plot of the normal displacement profile, (Equation |5]l, (d) surface 
plot of the electromagnetic energy functional, (Equation [6]l. In (d), there is significant optomechanical 
coupling density in the corner of the holes, where the crossbar meets the rail. Without the fillets, the field 
amplitude is concentrated in the corner and difficult to see. For this reason, the corners have been filleted to 
allow the optomechanical coupling density in the corners to be visualized. The fillets do not significantly 
affecting the optomechanical coupling (confirmed by simulation). 

The optomechanical coupling density can further be broken down into a mechanical part (the 
normal displacement profile) 

0m(r) = q-ft (5) 
and an optical part (the electromagnetic energy functional) 

0o(r)=A£|e|||^-A(£-i)|diP , (6) 

which can be separately visualized on the surface. This provides a quantitative method of assessing 
the separate optical and mechanical contributions and allows an intuitive approach to individually 
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engineering the optical and mechanical properties of the structure to enhance the optomechanical 
coupling of specific modes. 

Figurejij^a) shows the contribution to the optomechanical coupling, Lq^, of the breathing mode 
and fundamental optical mode from each "unit cell" of the structure. Summing the contributions 
from each unit cell yields LqI^. Figs. 4;b)-(d) show ^OM, 0m, and 0o, plotted on the surface of 
the nanobeam OMC for the fundamental breathing mode and the fundamental optical mode. In 
Fig.^b), it can be seen that there are two dominant and opposite contributions to the optomechan- 
ical coupling: one from the outside face of the rails and one from the inside face of the rails (in 
the comers of holes). Minimizing the cancellation between these two contributions is critical to 
achieving a small Lqm for the breathing mode (i.e. strong optomechanical coupling). The geom- 
etry of the "nominal" structure optimizes the coupling between the fundamental optical mode and 
the breathing mode, as shown below. 
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FIG. 5: For the fundamental breathing mode and the fundamental optical mode, (a) the dependence of 
the optomechanical coupling on the rail thickness (with oscillations in the data arising from accidental 
degeneracies with the cantilever modes), (b) the optical and mechanical mode profiles for rail thicknesses 
of 100 nm, 190 nm and 400 nm circled in red, green and blue respectively in (a), (c) comparison of the 
mechanical mode profiles when coupled (orange) and not coupled (purple) to cantilever modes, with the 
corresponding effect on in Lqm highlighted in (a). 



Since the breathing mode is drawn from a band edge at the F point, adjacent unit cells are 
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mechanically in-phase with each other and add constructively to the optomechanical coupling. 
This is in contrast to defect modes drawn from band edges at the X point, such as the pinch 
mode, where adjacent unit cells are mechanically out-of-phase, resulting in neighboring unit-cell 
contributions that tend to cancel. This cancellation reduces the optomechanical coupling unless 
it is specifically mitigated with extremely tight modal envelopes (see description of pinch mode 
optomechanical coupling below). 

The degree to which the different faces of the rails cancel each other's contribution to Lq^ is 
set by the attenuation of the optical field between the two edges, as the mechanical displacement 
of the two rails is fairly uniform. Thus, one would expect that varying the rail thickness, which 
changes the relative amplitude of the optical field on the two rail faces, would have a significant 
impact on the coupling. Figure [Sj^a) shows Lqm as a function of rail thickness, with Lqm of the 
nominal structure ( 1 90 nm rail thickness) circled in green. For rail thicknesses smaller than 1 90 nm 
(such as the 100 nm rail width, circled in red in Fig. [5]^a) and shown in Fig. [5]^b)), the amplitude 
of the optical field on the inside and outside edge of the field is becoming more and more similar. 
This results in a larger cancellation between the contributions to Lq^ on the inside and outside 
of the rails, decreasing the optomechanical coupling. This reasoning might lead one to believe 
that increasing the rail thickness should monotonically decrease Lqm (increase optomechanical 
coupling). However, for rail thicknesses larger than 190 nm (such as the 400 nm rail width, circled 
in blue in Fig. [5]^a) and shown in Fig. [Sj^b)), there is significant decrease in confinement of the 
optical mode because the light can partially "spill around" the holes through the wide rails. The 
mechanical mode, in contrast, stays relatively confined. The net effect is that the optical energy 
is "wasted" on parts of the structure that do not have significant motion, and the optomechanical 
coupling is again decreased. 

Just as Qm is affected by hybridization of the breathing mode with propagating and body modes, 
Lqm is affected by hybridization as q(r) is modified by the coupling to waveuide or body modes. 
This is responsible for the oscillations in Lqm seen in Fig. [Sj^a). The impact of coupling to the 
nanobeam body modes can be clearly seen in Fig. [Sj^c), where the breathing mode in a structure 
with a rail thickness of 230 nm has been plotted for two different beam lengths (number of total 
holes). For 47 total holes (circled in orange in Fig. |5^)), the breathing mode shape is altered 
significantly by the hybridization, causing the Lqm to deviate from the trend indicated by the red 
line in Fig. |5^). Shortening the structure by 2 holes (one on each side) decreases the coupling of 
the breathing mode to the propagating mode, returning Lqm to the trend line. 
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FIG. 6: For the fundamental pinch mode and the fundamental optical mode in the nominal structure, (a) 
FEM simulation of individual unit cell contributions to the total optomechanical coupling, (b) surface plot 
of the optomechanical coupling density, (c) surface plot of the normal displacement profile (Equation [5]), 
(d) surface plot of the electromagnetic energy functional (Equation|6]). 

The pinch mode is a localized, in-plane differential acoustic vibration. Each neighboring cross- 
bar vibrates 180 degrees out of phase with its nearest neighbors, since the pinch mode is drawn 
from a band edge at the X point. So although the optomechanical coupling contribution from 
each half of the structure (with respect to the y-z plane) is equal, such that the two halves add 
constructively to Lq^^, on either side of the y-z plane, contributions to LqI^ from neighboring 
crossbars tend to cancel. This puts a premium on mechanical localization, as a more localized 
pinch mode has a larger difference (and thus a reduced cancellation) between neighboring cross- 
bars. Although the envelope of the pinch mode's displacement profile is gaussian, each crossbar 
is very rigid, so the displacement of the compression and tension faces of each beam is essen- 
tially identical (but opposite). The gaussian envelope only serves to change the relative vibration 
amplitudes of neighboring crossbars. The optomechanical coupling contribution from each beam 
would be approximately zero if it weren't for the rapid variation of the optical mode's envelope, 
and the contribution of each crossbar to Lq^^ depends primarily on the difference in the optical 
energy density across the width of the beam. This would then lead one to believe that tighter 
localization, both optically and mechanically, would produce better optomechanical coupling for 
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FIG. 7: For the accordion mode and the fundamental optical mode in the nominal structure, (a) FEM 
simulation of individual unit cell contributions to the total optomechanical coupling, (b) surface plot of 
the optomechanical coupling density, (c) surface plot of the normal displacement profile (Equation [5]), (d) 
surface plot of the electromagnetic energy functional (Equation|6]l. 

this structure. Indeed, although the Lqm of the pinch mode in the structure shown is quite modest 
41 fim), Lqm can be reduced to less than 3 /urn by more tightly confining the optical and me- 
chanical modes by reducing the number of holes involved in the defect region. There is, however, 
a loss of optical Q associated with the increased confinement due to the larger optical momentum 
components associated with tighter spatial localization. However, the structure as shown has a 
radiation-limited optical Q greater than 10 million; so it can be quite reasonable to trade optical Q 
for higher optomechanical coupling. 

The last type of mechanical mode to be considered is the accordion mode (Fig. [7]). The rela- 
tively poor Lqm for the accordion mode in the nominal structure is partly due to the fact that the 
rails recoil against the motion of the cross bar, producing opposing optomechanical contributions 
within each unit cell. In addition, the coupling of the broad first order Hermite-Gauss envelope of 
the mechanical mode with the narrower optical mode induces cancellations in the optomechanical 
coupling contributions at the inflection points of the mechanical mode envelope. 

As discussed above, the accordion mode has a large effective mechanical bandgap. The dra- 
matically increased that results makes it worthwhile to investigate whether the structure can 
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FIG. 8: For the accordion mode with the fundamental optical mode, (a), the effective length as a function of 
total beam width, (b), individual unit cell contributions to the total optomechanical coupling for a structure 
with a beam width of 700 nm (circled in (a)), mode frequency of 3.97 GHz and effective motional mass of 
334 fg, with accompanying mechanical mode plot. The narrower mechanical mode (represented here by 
the deformation of the structure with color indicating relative strain) envelope results in drastically different 
optomechanical coupling contributions compared to Fig. |7] 



be modified to produce smaller Lqm- By reducing the width of the nanobeam, it can be seen from 
Fig. [8]^a) that the coupling is dramatically improved by almost two orders of magnitude when the 
width of the structure is reduced. As shown in Fig.[8]^b), for a beam width of 700 nm, the contribu- 
tions within each unit cell no longer cancel, due to the comparatively narrower mechanical mode 
envelope, and the structure yields Lqm = 3.67 ^um. In addition, simulations of the Qm show that 
the effective bandgap for this narrower structure is approximately 2 GHz, yielding an extremely 
large for a given number of holes in the mirror section (Q^ ^10^ for 35 total holes). It should 
be noted that the frequency of the accordion mode of the narrower structure is approximately 4 
GHz, up from 1 .5 GHz in the wider structure. 
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FIG. 9: (a) Fundamental optical mode of the double-heterostructrure OMC (geometry identical to that de- 
scribed in Ref.El), withXo « 1.5 i"in, 2rad ~ 2.7 X 10^, and Veff = 1.2 (Xo/w)^. (b) Breathing mechanical 
mode of the double-heterostructure OMC, with Vm = 9.3 GHz, and meff = 322 femtograms. (c) Optome- 
chanical coupling integrand plotted on the double-heterostructure OMC system's surface; the structure has 
an LoM = 1-75 jum for the optical-mechanical mode-pair from|9ja) and[9jb). 



VI. OPTOMECHANICAL COUPLING IN TWO-DIMENSIONAL OPTOMECHANICAL CRYS- 
TALS 

As a final example of how these methods can be used to understand the optomechanical cou- 
pling in periodic structures with complex mechanical and electric field profiles, we model a double 
heterostructure hexagonal photonic crystal slab resonator. This is a well-known optical system, 
which has been found to have radiation-limited quality factors in excess of twenty million, with 
experimental demonstrations exceeding quality factors of two millioiP^. The system consists of 
a hexagonal lattice of air holes in a silicon slab, with a single row of holes removed to create 
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a waveguide mode within the optical stop band (the defect pulls the waveguide mode from the 
conduction band); in addition, the spacing in the direction of the waveguide is abruptly decreased 
twice to provide longitudinal confinement. This structure is essentially equivalent to the nanobeam 
structure, with the optical and mechanical modes guided by Bragg reflection in the lateral direc- 
tion, as opposed to total internal reflection and hard boundaries in the nanobeam. With this in 
mind, we expect very similar optical and mechanical modes wherever lateral propagation out of 
the waveguide is prohibited by Bragg reflection. 

The fundamental optical cavity mode of the structure (the geometry is identical to that described 
in Ref.f^SI) has been reproduced by FEM simulation and shown in Fig. ^a). The structure also 
exhibits a lateral mechanical breathing mode at 9.3 GHz with a motional mass of 300 femtograms, 
modulating the width of the waveguide in a way that is analogous to the mechanical breathing 
mode of the nanobeam. The breathing mode displacement profile is shown in Fig.|9];b). 

Figure |9]^c) shows the integrand of the optomechanical coupling integrand (Eq. [3]) between 
the the optical mode and the mechanical breathing mode plotted on the surface of the structure. 
The structure is shown slightly tilted to allow the insides of the holes to be seen, which give the 
dominant contributions to the optomechanical coupling. It is interesting to note that the coupling 
comes almost entirely from the movement of a small part of the interior of the holes (i.e., the 
region of the inner sidewall of the hole, closest to the center defect region); this can be seen by 
comparing the top half of the structure to the bottom half (since the integrand is symmetric about 
the x-z plane). Since each row of holes provides an opposite contribution to its neighbors, it is 
necessary to have a rapidly decaying optical envelope to achieve small Lqm, which is the case for 
the optical mode shown here. There is also a very small, opposing contribution from the center 
waveguide due to buckling/extrusion of the structure as the width is modulated. Just as in the case 
of the nanobeam, this optical-mechanical mode -pair has a very strong dispersive coupling, and 
evaluating the integral yields an effective length of only 1.75 /jm. 

Optically, the structure has a complete photonic bandgap for in-plane propagation, but, with 
a hole size to lattice constant ratio of r/A = 0.26, there is no corresponding in-plane mechanical 
bandgap. This makes the structure susceptible to mechanical loss mechanisms similar to those 
of the nanobeam. However, the two-dimensional hexagonal lattice, as well as other two dimen- 
sional Bravais lattices, can have simultaneous optical and mechanical bandgapd^^'^, allowing the 
possibility of highly localized, low-loss optical-mechanical mode -pairs with very small effective 
lengths and motional masses. 
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